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ABSTRACT 


Employing  a  recent  homogeneous  strain-dependent  creep  damage 
theory,  the  propagation  of  a  failure  front  in  a  beam  under  pure  bend¬ 
ing  is  studied.  A  local  inhomogeneous  strain-dependent  creep  damage 
theory  is  then  postulated,  based  on  creep  damage  data  obtained  at 
J.R.C.  Ispra  for  specimens  of  various  lengths.  Using  this  inhomo¬ 
geneous  theory,  solutions  for  various  loading  conditions  are  obtained, 
and  are  shown  to  be  in  better  agreement  with  observation  than  previous 
theories.  Finally,  the  effect  of  random  material  parameters  on  the 
rupture  time  is  considered  for  the  constant  tensile  stress  test  and 
for  the  beam  under  pure  bending. 
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1.  INTRODUCTION 


With  the  development  of  heat  resistant  metallic  alloys,  it  is  now 
possible  to  increase  the  operating  temperatures  in  thermal  machines 
(e.g.,  jet  turbines  and  nuclear  reactors)  and  thereby  obtain  better 
thermodynamic  efficiencies.  At  these  elevated  operating  temperatures 
the  phenomenon  of  creep  is  often  very  significant,  and  in  particular 
the  problem  of  creep  rupture  is  of  major  concern  when  considering  the 
safety  of  such  thermal  machines.  Creep  rupture  in  even  a  small  struc¬ 
tural  component  of  a  nuclear  reactor  may  not  only  jeopardize  the 
operation  of  the  entire  system,  but  can  also  result  in  a  serious  hazard 
to  the  surroundings.  With  the  availability  of  high  speed  computers 
and  modern  numerical  techniques  the  designer  may  now  obtain  accurate 
creep  rupture  analyses ;  this  accuracy  is  primarily  limited  by  the 
reliability  of  the  creep  rupture  theory  employed.  The  present  research 
is  directed  toward  the  development  of  improved  creep  rupture  theories, 
through  the  inclusion  of  the  effects  of  strain,  inhomogeneity  and  ran¬ 
domness,  and  the  application  of  these  theories  to  illustrative  examples 
Deterministic,  homogeneous,  uniaxial  creep  rupture  theory  has  been 
developed  along  two  major  approaches:  I.  The  continuum  mechanical 
approach,  and  II.  The  physical  metallurgical  approach.  The  first 
approach  has  been  used  by  researchers  specializing  in  continuum  and 
structural  mechanics,  and  in  applied  mathematics.  For  the  most  part 
such  researchers  have  been  interested  in  the  macroscopic  phenomenon  of 
rupture  due  to  creep,  and  have  sought  an  interpretation  of  the  charac¬ 
teristics  of  uniaxial  creep  strain  on  rupture  data  curves.  The  second 
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approach  has  been  used  by  researchers  specializing  in  metallurgy  and 
solid  state  physics.  They  have  been  interested  in  quantitative  and 
qualitative  studies  of  the  microscopic  damage  which  leads  to  rupture. 
Multiaxial  creep  rupture  theory  is  then  usually  obtained  by  a  ten- 
sorially  reasonable  extension  of  uniaxial  theory  with  the  damage 
itself  usually  still  treated  as  a  scalar,  although  in  some  more 
general  formulations  the  damage  lias  been  treated  as  a  tensor. 

In  the  present  study  we  first  discuss  various  deterministic, 
homogeneous,  strain-independent,  uniaxial  and  multiaxial  creep  rupture 
theories.  Then  the  deterministic,  homogeneous,  uniaxial  tensile, 
strain-dependent  creep  damage  theory  recently  developed  by  Belloni, 
Bemasconi,  Cozzarelli  and  Piatti  [1,2,3],  using  both  the  continuum 
mechanical  and  physical  metallurgical  approaches,  is  discussed  in 
detail,  and  generalized  to  the  cases  of  compressive  and  multiaxial 
stress.  The  problem  of  the  development  of  a  failure  front  in  a  beam 
under  pure  bending  is  studieo  using  this  strain-dependent  creep 
damage  theory.  A  local  damage  law  is  then  obtained  for  inhomogeneous 
materials  under  constant  tensile  stress,  by  means  of  an  analysis  of 
creep  damage  data  for  specimens  of  various  lengths,  and  this  law  is 
extended  to  the  case  of  variable  tensile  or  compressive  stress.  This 
local  damage  law  for  variable  s*  '■ess  is  used  to  obtain  solutions  in  a 
bar  under  constant  load  accompanied  by  lateral  contraction.  Finally, 
for  the  problems  previously  treated  deterministically,  the  effect  of 
randomness  in  the  material  parameters  on  the  rupture  time  is  analyzed 
by  probability  techniques. 

Section  2  contains  the  literature  survey  of  previous  creep 


rupture  theories.  The  strain-dependent  creep  damage  theory  and  its 


application  to  a  beam  under  pure  bending  is  presented  in  Section  3. 

In  Section  4  we  obtain  the  local  damage  law  for  inhomogeneous  materials, 
and  solutions  are  obtained  for  various  loading  conditions.  Finally, 
the  randomness  of  the  material  parameters  is  considered  in  Section  5, 


and  an  overall  summary  is  given  in  Section  6. 
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2.  LITERATURE  SURVEY  OF  PREVIOUS 
CREEP  RUPTURE  THEORIES 

As  mentioned  before,  there  are  two  major  approaches  to  the  develop¬ 
ment  of  uniaxial  creep  rupture  theories,  i.e.,  continuum  mechanical 
and  physical  metallurgical.  An  early  accomplishment  in  the  continuum 
mechanical  approach  was  due  to  Hoff  [4],  who  assumed  a  ductile  mode  of 
creep  rupture.  This  phenomenon  occurs  at  high  tensile  stresses  as  the 
cross-sectional  area  A  narrows  in  time  due  to  lateral  contraction  and 
theoretically  decreases  to  zero  at  rupture;  such  ductile  rupture 
generally  occurs  in  a  short  period  of  time.  The  logarithmic  strain- 
displacement  relation  for  large  deformation,  along  with  an  assumption 
of  incompressible  lateral  contraction  for  a  bar  under  constant  load, 
yields  the  relation 


c  =  o  e 
o 


(1) 


where  o  is  the  time-dependent  stress,  oq  is  the  initial  stress, 
and  c ^  is  the  time-dependent  creep  strain.  In  addition,  Hoff 
employed  the  Norton's  law  for  steady  creep 


£  =  Ko 

c 


(2) 


where  K  and  n  are  material  creep  parameters.  Employing  the  ductile 

rupture  criterion  (A  -*■  0  and  t..us  e^  -*■  00  at  rupture)  together  with 

Equations  (1,2),  Hoff  obtained  an  expression  for  the  ductile  rupture 

time  t  as 
R 
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t 


R 


1 


n  K  a 


n 


(3) 


o 

which  exhibits  a  linear  behavior  when  plotted  as  In  oq  versus  In  tR. 

In  practice,  materials  are  frequently  subjected  to  tensile  stresses 
below  that  required  for  ductile  rupture,  and  it  takes  a  longer  time  for 
a  bar  to  fail  and  it  does  so  without  obvious  lateral  contraction.  This 
mode  of  failure  is  called  brittle  creep  rupture.  In  order  to  describe 
this  type  of  rupture,  Kachanov  [5]  accounted  for  material  deterioration 
(or  damage)  by  introducing  a  quantity  called  the  continuity  tp  in 
accordance  with 


ip  = 


(4) 


where  Aq  is  the  initial  cross-sectional  area  and  Ar  is  the  effec¬ 
tive  undamaged  cross-sectional  area.  Since  the  damage  is  assumed  to 
be  distributed  over  the  material  rather  than  concentrated  at  a  crack, 
this  approach  is  sometimes  called  a  distributed  damage  theory  (see 
Krajcinovic  [6]).  Instead  of  the  continuity  it  is  generally  more 
convenient  to  employ  the  material  damage  parameter 


U)  =  1  -  ij>  = 


A  A 
o  -  r 


(3) 


and  Kachancv  postulated  for  the  damage  rate  the  power  law 


u  =  B 


o 

1  -  u 


(6) 


where  B  and  v  are  material  damage  constants.  If  we  ignore  lateral 

contraction  under  constant  load,  o  is  simply  equal  to  the  constant 
P 

stress  a  =  —  where  P  is  the  constant  load.  Equation  (5)  indicates 

o  A.  o 
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Chat  if  we  assume  that  there  is  no  damage  in  the  initial  state  then  at 
t  =  0  ui  =  0  (or  =  A0)  ,  and  also  that  at  rupture  w  -*■  1  (i.e., 

Aj.  =  0)  although  this  theoretical  limit  is  never  actually  reached. 
Using  this  brittle  rupture  criterion  and  Equation  (6) ,  the  rupture 
time  for  purely  brittle  creep  rupture  at  constant  stress,  is  obtained 
as 


R  B (IK)) o  V 
o 

This  relation  again  shows  a  linear  behavior  in  a  In  oQ  versus  In  t^ 
plot . 

Kachanov  also  used  Equations  (1)  and  (2)  for  a(t)  in  Equation 
fb)  for  the  general  constant  load  test  with  lateral  contraction  included. 
This  is  called  a  one-way  damage  creep  interaction,  since  deformation 
is  assumed  to  affect  the  rate  of  damage,  while  damage  is  assumed  not  to 
affect  the  rate  of  deformation.  The  rupture  time  is  obtained  from 
this  approach  as 


n 


(8) 


where  t  and  t  are  the  Hoff  and  Kachanov  rupture  times  given 
Kli  KK. 

bv  Equations  (3,7)  respectively. 

Following  Hoff's  and  Kachanov's  basic  ideas,  many  researchers 
have  introduced  improvements  in  either  the  creep  law  and/or  the 
damage  law.  Rabatnov  [7]  proposed  the  two-way  damage  creep  interac¬ 


tion  equations 
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K  on(l  -  w)_q 

(9) 

B  oV(l  -  u)'r 

(10) 

where  q  and  r  are  material  powers,  and  where  the  presence  of  w 
in  Equation  (9)  may  be  interpreted  as  an  inclusion  of  tertiary  creep 
since  increase  with  time.  From  these  equations,  one  now  obtains 
the  rupture  time  simply  as 


CR  = 


B(l-*r)o 


(11) 


and  in  addition  the  creep  strain  at  rupture  is  obtained  as 


1  K  n-v 
eR“l+r-qB0o 


(12) 


Odqvist  [8]  added  a  term  to  the  creep  law  to  account  for  nonlinear 
instantaneous  deformation  (and  for  transient  creep  in  an  approximate 
manner),  i.e., 


e 

c 


=  Ko 


)  +  K  o 


(13) 


where  Kq  and  nQ  are  instantaneous  deformation  material  parameters. 
Broberg  [9]  also  added  a  term  to  Rabotnov's  Equation  (10)  to  account 
for  instantaneous  damage,  and  proposed  the  modified  pair  of  equations 


=  G'(S)S  +  F (S) 

(14) 

-  g’(S)S  +  f(S) 

(13) 

where  S  ■  a/(l  -  w) .  BostrOn  et  al.  [10]  proposed  a  logarithmic 
definition  of  damage  by  analogy  with  the  logarithmic  strain  as 
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A 

*  -  In  (16) 

A 

r 

Westlund  [ 1 1 J  and  Belloni  et  ai.  [12]  have  given  interesting  compari¬ 
son  of  the  various  theories  developed  from  the  Hoff  and  Kachanov 
approaches.  In  the  next  section  we  shall  present  in  detail  the  recent 
strain-dependent  creep  damage  theory  of  Belloni,  Bemasconi,  Cozzarelli 
and  Piatt i  [  1,2,J] . 

In  addition  to  the  developments  following  Hoff's  and  Kachanov's 
damage  theories,  there  are  aiso  ■ ome  approaches  from  the  energy  point 
of  view.  Kopecki  and  Walczak  [13]  used  the  dissipated  energy  as  the 
rupture  criterion,  i.e.,  at  rupture 

f(EQ)  =  constant  =  K  (17) 

where  is  the  specific  dissipated  energy  and  K  is  a  material 

rupture  constant.  It  is  postulated  in  [13j  that  if  the  energy  barrier 
as  expressed  by  Equation  'l-'1  has  been  overcome,  the  material  loses 
its  ability  to  dissipat.  energy  and  rupture  occurs.  Janson  and  Hult  [14 
employed  a  combined  approach  wit.  a  single  defect  along  with  distributed 
damage,  and  the  strain  energy  was  expressed  in  terms  oi  an  effective 
stress  S  .  it  is  postulated  in  [14]  that  when  this  energy  reaches  a 
critical  value  equal  to  the  work  needed  to  create  two  new  separate 
surfaces  (i.e.,  a  single  crack)  the  material  will  rupture. 

The  physical  metallurgical  approach  has  been  employed  by  Hult  et  al 
[15],  where  an  attempt  was  made  to  relate  damage  to  the  change  in  the 
velocity  of  sound  in  the  material.  Also,  the  variation  of  electrical 
resistivity  was  used  as  an  index  Jr  ap  .  Furthermore,  neutron 
diffraction  or  x-ray  diffract  on  techniques  have  also  been  used  to 
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measure  damage.  Piacci  et  al.  [lt>]  Lave  developed  very  sensitive 
differential  density  measurement  techniques  as  an  index  of  damage  due 
to  void  nucleation  and  growth,  and  this  is  discussed  further  in 
Section  3.  Finally,  the  tedious  method  of  actual  void  counting  with 
the  aid  of  optical  and  electron  microscopes  has  also  been  employed. 

Due  to  the  complexity  of  creep  rupture  under  multiaxial  states  of 
stress,  there  is  at  present  no  completely  satisfactory  three- 
dimensional  theory.  In  fact,  even  the  elementary  one-dimensional 
theory  is  complicated  by  the  fact  that  much  more  damage  occurs  under 
tension  than  under  compression  (e.g. ,  see  Schiller  et  al.  [17]). 
Odqvist  [18]  generalized  the  uniaxial  theory  by  replacing  the  uniaxial 
stress  in  the  damage  law  by  the  maximum  principal  tensile  stress  in 
accordance  with 


where  o  is  the  maximum  principal  tensile  stress,  and  t,,  is  the 

max  r  h 

rupture  time  under  the  constant  uniaxial  stress  o  .  However, 

K 

Johnson,  Henderson  and  Khan  [IV]  and  Henderson  and  Mathur  [20,21]  have 
experimentally  demonstrated  tiiat  the  rupture  time  of  different 
materials  appears  to  depend  on  different  rupture  criteria,  liayhurst  [22] 
has  shown  that  copper  and  aluminium  alloys  demonstrate  two  extreme  types 
of  rupture  behavior  under  multiaxial  stress.  For  copper  a  maximum  ten¬ 
sile  stress  criterion  such  as  i.quation  (18)  was  found  to  be  valid, 
while  for  aluminium  alloys  a  maximum  shear  stress  criterion  was  more 


appropriate. 
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Martin  and  Leckie  L  2 j)  ]  and  nnyhurst  and  Leckie  [24]  have  proposea 
a  combined  theory  which,  included  these  two  limiting  cases.  For 
homogeneous  damage  (i.e.,  t lie  eitect  of  damage  is  included  in  all  three 
principal  stresses)  a  potential  lunction  $  was  assumed  as 


} n+1  ! o  -  o 

o  1  J  i  l 


(19) 


ihe  principal  strain  rates  were  derived  from  (  as 


Hi 


1  1 


(20a) 


t 


I  i 


u 


(20b) 


where 


il 


ana 


li 

t. 

0 


tII 

u 

o 


and  t 
are  norr. 


Ill 

inal 


are  the 
values . 


,  I  >  )  are  tke  principal  stresses, 

cor  responding  principal  strain  rates,  and 
Thou  the  damage  law  was  given  by 


j 


(21) 


and  it  follows  that  for  o  =  j  Liquation  (21)  corresponds  to  a  maximum 
tensile  stress  criterion  and  for  u  **  1  Liquation  (21)  corresponds  to 
a  maximum  shear  stress  criteii  'or  noni.omogeneous  damage  (i.e., 

the  effect  oi  damage  is  included  only  in  o  )  the  potential  function 
was  given  b> 


=  l 
(it 

^  o  ' 


1_ 

I  n+1 


1* 


ilij 


(22) 


The  principal  strain  rates  then  were  derived  from  Lquation  (22)  as 


11 


-eTTT  =  e 
III  o 


III 


/o. 


(23a) 


e 


II 


=  0 


(23b) 


and  the  damage  law  was 


°I 

*  =  -b  ~  -  e  o 
Ul* 


hi 


“I  v 


/a. 


(24) 


The  previous  multiaxial  theories  express  the  material  damage 
in  terms  of  scalar  quantities.  On  the  other  hand,  Murakami  and  Ohmo 
[25]  have  assumed  that  material  damage  can  be  expressed  as  a  symmetric 
tensor  of  rank  two.  Recently,  a  tensor  of  rank  eight  was  proposed  by 
Chaboche  [26]  for  the  multiaxial  material  damage.  This  damage  tensor 
was  reduced  to  a  tensor  of  rank  two  using  the  symmetry  of  the  stress 
tensor.  In  the  following  section  we  present  a  recent  damage  theory, 
in  which  strain  plays  a  significant  role. 
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3.  THE  STRAIN-DEPENDENT  CREEP  DAMAGE  THEORY 

In  cue  previous  section  all  of  the  damage  laws  were  expressed 
as  independent  of  the  strain,  but  recent  experimental  results  obtained 
by  belloni,  Bernasconi  and  Piatt i  at  the  Euraton  Laboratory  [lj  sug¬ 
gest  that  strain  may  in  lact  be  a  dominant  variable  in  damage  accumula¬ 
tion.  These  results  establish  the  connection  between  relative  density 
variation,  which  is  taken  as  a  measure  of  damage,  and  the  creep  strain, 
temperature,  stress  and  time.  Thus,  the  damage  law  is  in  general 
expressed  as 

— =  f  U  ,T,o,t)  (25) 

P  c 

o 

where  o  is  tiie  initial  density,  -do  is  the  reduction  in  density 
due  to  void  ormation  and  growth,  and  T  is  the  temperature . 

j.a  Strain-Dependent  creep  Damage  Power  Law 

The  strain-dependent  creep  carnage  theory  presented  in  [lj  has  been 
experimental J v  verified  using  verv  sensitive  dilieiential  density 
measurement  techniques  | 27-dVj  .  By  these  techniques  the  damaged  speci¬ 
men  and  an  undamaged  dummy  ar  >.  d  Li  t  n  in  air  and  in  a  heavy 

liquid  using  an  analytical  balance,  and  then  relative  density  variations 
as  small  as  10  J  or  even  lo  ’  are  computed  from  a  simple  formula. 
by  using  an  extensive  conipu; rr  tat i«t leal  regression  analysis,  a  good 
agreement  was  established  witn  l  he  experimental  i  esolts  for  the  case 
o;  one  dimensional  constant  ten  cress  tas  approximated  by  constant 

loan),  using,  the  separable  ;  .wei  law  luim  tor  the  damage  law 
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Ap  _  a  -3/T  Y  6 

- -Ce  e  at 

p  o  c  o 


(26) 


Equation  (26)  simplifies  to 


Ap  _  a  Y  <5 

-  — 1 -  =  c  e  a  t 

P  co 

o 


(27) 


if  the  temperature  is  constant  at  Tq  and  C  is  the  constant 

C  =  CQe  o  .  In  Equations  (26,27)  a,  0,  Y,  6  are  exponents  which 

characterize  respectively  the  effects  of  creep  strain,  temperature, 

stress  and  time  on  the  damage.  The  experimental  data  indicate  that 

the  creep  strain  is  always  the  most  important  factor  among  these  four 

damage  variables.  Also,  the  four  exponents  were  shown  to  be  positive 

constants,  which  depend  on  the  material  and,  within  certain  limits,  do 

not  depend  on  the  test  conditions.  Also,  in  [1]  the  Kachanov  damage 

parameter  w  was  shown  to  be  essentially  equal  to  -  by  means  of 

Po 

a  simple  conservation  of  mass  calculation. 

In  the  present  study,  Equation  (27)  will  be  rewritten  in  the  more 
general  form 


D(t)  =  C  e  a  o  Yt^ 
c  o 


(28) 


where  D(t)  denotes  the  damage  and  does  not  necessarily  correspond 

with  .  Combining  Equations  (2,28)  (i.e.,  excluding  transient 

Po 

creep)  one  obtains  the  rupture  time  for  constant  tensile  stress  as 


tR 


(C  Ka0ona+\ 


1 

a+6 


(29) 


where  D  is  the  critical  value  of  damage  at  rupture.  This  is  again 
K 

a  linear  relation  in  a  In  oo  versus  In  tR  plot.  By  assuming  that 


d 


u 


the  rupture  times  in  Equations  are  the-  sane,  one  car.  estu  -is-, 

a  correspondence  between  Kachanov’s  theory  and  Equation  1.28),  yielding 


na+\ 

u+6 


B 


_1_ ,  C_  K_ ,  ' 

1+V>K- 


OD 


This  is  an  important  result  since  it  makes  it  possible  to  interpret 
Kachanov's  theory  from  the  metallurgical  point  of  view. 

Cozzarelli  and  Bernasconi  12]  extended  Equation  (28)  to  the  case 


of  variable  tensile  stress  c(t)  as 

t 

D(t)  -  C  ec(t)a  ;  c(t'/,/6dt’"i 


(32) 


o 

where  for  o(t)  =  Equation  t32)  reverts  to  the  form  in  Equation  (28). 
Furthermore,  it  possesses  the  physically  necessary  characteristic  that 
D(t)  be  continuous  even  if  the  stress  history  is  discontinuous.  Using 
Equation  (32),  rupture  times  v«.-re  also  determined  for  various  loading 
conditions,  such  as  step-stress,  constant  stress  rate,  relaxation  and 
constant  load  with  lateral  cont :  action.  Hult  [30]  and  Chrr. anovski  [31] 
assumed  that  there  is  no  dar.age  produced  under  compression,  and  in  [17] 
it  was  shown  that  the  experimental  data  essentialiv  justify  such  an 
assumption.  Hence,  Equation  a  oc  peneralized  further  as 


i 

'j(t)  =  C  i  (t)u-  >  (t  ’  '  :  l[o(t  *)]dt  ’  l  (33) 

c  U  J 

u 

where  o(t)  can  be  either  tensile  or  eorprc-oSive ,  and  U(o)  is  the 


unit  step  function. 


under  pure  bending. 
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J.’d  Damage  in  a  Beam  under  Pur*  Bending 

Consider  an  initially  straight  beam,  with  a  uniform  rectangular 
cross-section  of  width  w  and  initial  effective  depth  2hQ ,  and  sub¬ 
jected  to  pure  benuing  by  constant  moment  M  .  As  mentioned  in  the 
last  section,  there  is  essentially  no  damage  produced  under  compression 
and  therefore  the  beam  will  deteriorate  only  on  that  side  of  the 
neutral  axis  that  is  under  tensile  stress.  Since  the  maximum  tensile 
stress  occurs  at  the  top  outside  fiber  (.see  Figure  1),  rupture  will 
initiate  at  this  fiber  and  then  a  failure  front  (.see  [5])  will  propagate 
down  through  the  entire  cross-section.  The  neutral  axis  is  at  the  cen¬ 
troid  of  the  undamaged  material,  and  we  denote  the  initial  neutral 

axis  as  x  and  the  initial  axis  in  the  direction  of  the  depth  as  v  . 

o  o 

As  rupture  develops  and  the  failure  front  propagates  the  neutral  axis 
moves  down  and  is  now  identified  as  x  ;  the  intersection  of  x  with 
y  defines  the  moving  origin  of  the  instantaneous  depth  coordinate  y  . 

In  order  for  Norte;’.'.,  law  1)  to  be  valid  for  both  tensile  and 
compressive  stresses,  it  is  rewritten  in  the  form 

t  =  k!o!  sgn(.o)  (37) 

c 

where  sgnfo)  is  the  signum  function  defined  as 


i 


S(..ihc)  -  (3b) 

iy-1  lor  o  •  u 

and  wht-ru  n  is  not  necessar  il>  -••str  i  cted  to  ole  posit  -e  integers 
Kmploving  the  Fuler-Bernoul  1  i  a.  eumpt  if .>!  plan*.-  .  ections  remaining 
plane,  we  obtain  the  strain  rat*  in  terms  of  the  curvatui e  t  at  * 

as 


An  expression  for  the  bending  moment  is  then  found  from  Equations 
(37,39)  as 


M  = 


oydA  =  l]“|n  3  sgn(x) 


(40) 


where 


J 


dA 


(41) 


which  depends  on  the  geometry  of  the  cross-section  and  the  creep 
power  n  .  For  a  beam  of  rectangular  cross-section  with  width  w  and 
instantaneous  depth  2 h  ,  we  have 


J 


-iS_  wh 
2n+l 


(42) 


2  3 

For  n  =  1  this  reduces  to  -wh  ,  which  is  the  ordinary  second 
moment  of  area  for  a  rectangular  beam.  The  stress  is  obtained  from  a 
combination  of  Equations  (37,39,40)  as 
1 

°  =  1  j~‘~ '  sgn(y)  (43) 


Initially,  the  maximum  tensile  stress  o  in  the  beam  occurs  at  the 

7  max 


top  outside  fiber  y  =  yQ  =  hQ  ,  and  is  given  by 


H(ho> 


max 


(44) 
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2+  - 

,  .  2n  n 

where  J  =  - - ~wn 

o  2n+l  o 

The  entire  cross-section  of  the  beam  will  resist  the  applied  bend¬ 
ing  moment  until  a  failure  front  forms  at  v  =  h  .  Since  the  maximum 

'  o 

tensile  stress  will  remain  constant  up  to  the  time  of  initial  rupture, 
this  time  is  given  by  Equation  f2i>)  as 


1 


,  n 


■R1 


f  D^  j-M(,h  )*-i  -fnu+yj;a+6 

j-  ; 

,t.  h  ^  o  -  } 


(45) 


Beginning  at  t  -  t  the  failure  front  will  travel  downward ,  and  the 
maximum  tensile  stress  is  then  at  the  failure  front  and  travels  along 
with  it.  During  this  rupture  propagation  process  the  maximum  tensile 
stress  is  no  longer  a  constant.  Figure  i  shows  that  at  any  time 
t  >  tR-  the  thickness  of  the  deteriorated  region  is  24U),  and  the 
effective  depth  of  the  beam  has  been  reduced  to  a  value  2h(t)  <  2h0  . 
It  is  clear  l'rcm  the  figure  that  hQ  =  lift)  +  SU)  and  y(t)  =  v  (t) 

+  h  -  ii(t)  ,  where  v(t)  is  ;  .usured  from  the  shifted  position  of 
the  neutral  axis.  ihe  nistribution  oi  stress  is  then  given  by 


1 
■  n 


.  /  ..  n  _  M  I  y  ( t )  : 

1  - ,  «gmy; 


J(t) 


1.40.1 


and  the  strain  is 


(t)  =  K  i  o  (t '  )  j  ‘  ::gn(.o)dt ' 


(47) 


From  Equations  (33,46,47),  ve  finally  obtain  the  damage  at  y(t;  >  0 


b(t) 


n  :  •  dt.v  jv/A  'ivtJV(t,)]V4,i 


f  }  [  V  ( t  )+h  —lift  '  )  | 

Cf  KM 


■Hf) 


n 


( 


dt 


(.48) 


19 


Also,  at  the  time  t  >  tD1  ,  with  the  failure  front  at  the  position 
y  =  h(t)  ,  it  is  clear  that  at  this  front  y^(t)  +  hQ  =  2h(t)  and 
D(t)  =  D  .  Therefore,  Equation  (48)  yields 

I\ 


D  =  C  K  M 
K 


'2n+l' 


f2n+l' 

n 

w' 

[2nw  j 

2nw 


V/6r 


[2h(t)-h(t')h(t')“1"2n  dt' 

_  X.  _  lx 

[2h(t)-h(t')]h(t')  6n  6  dt' 


(49) 


From  the  experimental  data  given  in  [3],  we  find  that  although  in 
general  y  ^  6n  the  value  of  y  is  in  the  vicinity  of  fin  .  For 
the  sake  of  simplicity  we  shall  assume  here  that  y  =  fin  (also  see  [2]), 
in  which  case  Equation  (49)  then  simplifies  to 


C  K  M 


a  ncrf  yl 


(J 

[  2n+l] 

n 

\ 

[2nw  j 

• 

[2h(t)-h(t')]h(t') 


-l-2n 


dt’ 


a+6 


(30) 


Differentiating  Equation  (50)  with  respect  to  t  ,  we  find  that 


2h 


h(t') 


-l-2n 


dt  '+h 


-2n 


=  0 


(31) 


where  h  designates  dh/dt  ,  In  addition,  we  have  the  initial  condi¬ 
tion 


h  =  h  at  t 
o 


■RI 


(52) 


Upon  substituting  Equation  (52)  back  into  Equation  (51)  ,  another 
initial  condition  is  found  as 

A  =  "  2l^  at  C  =  Si  (53) 
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A  second  order  differencial  equation  lor  h  n.av  be  obtained  by 
differentiating  Equation  (51)  with  respect  to  t  and  eliminating  the 
integral,  yielding 

h  +  2(n-l)  h  ^fh)“  =  U  (54) 


This  equation  is  of  the  same  form  as  that  derived  by  Kachanov  [5]  and 
Odqvist  [la].  Solving  differential  Equation  (54)  with  initial  condi¬ 
tions  (52,53;  ,  we  then  obtain  the  useful  result 


1  + 


2n-l 


1 


1  li_ 


ln-1- 


(55) 


At  the  instant  of  complete  rupture,  i.e.,  when  h  =  0  ,  the  rupture 
time  is  given  simply  as 


1  +  - 


2n-l 


(56) 


At  the  instant  that  the  ^oss-section  has  deteriorated  to  half 

of  its  initial  effective  depth,  i.e.,  h  =  yh^  ,  we  obtain  the  half 

rupture  time  t  ,  as 
K 


tR'  ,  ,  2  -2n+l, 

T~  '  1  +  2n-l  1 

KI 

V 


(57) 


The  ratio  as  computed  fro.  Lu  at  ion  (57)  is  only  slightly  dif- 

lk;  tR 

ferent  from  — -  [Equation  (50)],  and  this  difference  decreases  as 

Kl 

n  increases.  This  interesting  result  indicates  that  it  will  take 
only  a  small  part  of  the  total  :iretime  for  the  cross-section  of  the 
beam  to  go  from  a  state  of  502  deterio’  t  ion  to  complete  rupture,  and 


the  rupture  strength  of  the  team  decreases  as  n  increases.  Figure  2 


shows 


plotted  as  a  function  of 


for  the  following  two  steels 


[3]: 


AISI  310  stainless  steel  at  600°C,  where  n  =  7 
2,25Cr  IMo  ferritic  steel  at  550°C,  where  n  =  5.6 


The  figure  clearly  shows  the  acceleration  of  the  propagation  of  the 
failure  front  and  the  weakening  of  the  beam  as  n  increases.  As 
indicated  in  [18],  when  half  the  depth  of  the  beam  has  deteriorated, 
the  rupture  process  becomes  more  and  more  ductile  due  to  the  develop¬ 
ment  of  a  high  stress  distribution.  Since  the  above  development  does 
not  contain  a  Poisson  effect  such  as  in  the  ductile  rupture  theory  of 
Hoff,  it  is  less  accurate  in  this  latter  time  interval. 

In  the  next  section,  the  inhomogeneity  of  the  material  will  be 
considered  in  order  to  further  improve  the  strain-dependent  creep 
damage  theory. 
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A.  INHOMOGENEOUS  STRAIN-DEPENDENT  CREEP  RUPTURE 

The  theories  in  the  previous  two  sections  all  assume  that  the 
material  is  homogeneous ,  and  accordingly  the  material  parameters  are 
independent  of  position.  However,  in  actual  materials  there  are 
numerous  imperfections,  which  results  in  significant  inhomogeneity 
in  the  structure  sensitive  creep  parameters.  It  was  noted  by  Broberg 
[.33]  that  both  the  Hoff  and  Kachanov  theories,  as  applied  to  a  bar  in 
tension,  have  a  tendency  to  overestimate  the  lifetime  to  rupture.  He 
showed  that  this  may  be  attributed  to  variations  of  the  creep  strain 
and  creep  damage  parameters  (K  and  B)  along  the  axis  of  the  bar  result¬ 
ing  from  inhomogeneity  in  the  material.  Similarly,  we  shall  clearly 
demonstrate  in  this  section  that  in  the  strain-dependent  creep  damage 
theory  the  inhomogeneity  of  the  material  (coupled  with  the  fact  that 
the  experimental  values  of  the  critical  damage  D  are  measured  over 
a  finite  distance  rather  than  at  a  point)  can  account  for  the  fact  that 
D^  is  far  less  than  the  theoretical  limit  value  of  1. 

A. a  Inhomogeneous  Strain-Depenui nt  Damage  Law  foi  Constant  Stress 

The  problem  of  spatial  f  h.  'uns  (generally  random)  in  the 

material  creep  parameters  is  discussed  in  detail  in  Cozzarellt  [JA], 
where  it  is  noted  that  for  simplicity  one  may  in  some  cases  neglect 
spatial  tluctuation  in  the  creej  power  n  while  allowing  -nly  tin. 


creep  reciprocal  viscosity  K  to  fluctua'e  (as  i.  broberg  [33]  and 
Broberg  and  Westlund  [33]).  Ex  ten.  ,g  th-.s  concept  to  the  present 


23 


creep  rupture  problem,  we  assume  that  all  of  the  powers  n,  a,  y,  6 
are  constant  while  K  and  the  damage  coefficient  C  may  fluctuate 
spatially  due  to  material  inhomogeneity.  For  one  dimensional  constant 
tensile  stress,  the  Norton  creep  law  then  takes  the  form 

ec(x,t)  =  K(x)oont  (38) 

where  K(x)  depends  on  the  axial  coordinate  x  .  In  order  to  account 
for  inhomogeneous  damage  that  is  present  prior  to  loading,  we  intro¬ 
duce  tne  concept  of  a  damage  difference  D  ,  which  is  defined  as 

D(x,t)  =  DT(x,t)  -  Dq(x)  (59) 

where  DT(x,t)  is  the  total  damage  and  Dq(x)  is  the  initial  damage 
at  t  =  0  .  At  a  particular  constant  temperature  Tq  ,  we  postulate 
by  extension  of  strain-dependent  law  (28)  the  local  damage  law 

D(x,t)  =  C(x)ec(x,t)0‘ooYt'5  (60) 

Upon  substituting  Equation  (58)  into  Equation  (60)  ,  we  then  obtain 

Dfa.O  -  HWo^-V-4  (61) 

where  we  have  introduced  the  combined  material  parameter  H(x)  = 
C(x)K(x)a. 

For  a  particular  applied  constant  tensile  stress  o at  constant 

temperature  Tq  ,  a  one-dimensional  specimen  of  length  l  will  reach 

rupture  at  some  position  *  and  at  a  particular  time  t  .  In  the 

O  OK 

neighborhood  of  this  position  we  may  compute  the  average  damage  at 
-  AV 

rupture  D  (x)  over  interval  *  i  x'  $  *  as 

K  O 


X 


j  VX,’tCR;dX' 

D  AV)  =  A— - 

R  '  '  X  -  X 


ll(x')dx' 


x 

o 


X  -  X 

o 


Q 


2A 


(62) 


where  Q  =  0  na+^t  a+^'  #  The  experimental  results  given  in  [17] 
contain  very  useful  data  for  the  fractional  density  change  plotted  as 
a  function  of  the  mass  for  various  one  dimensional  specimens  of  dif¬ 
ferent  lengths  cut  from  the  test  specimen  at  the  position  of  rupture. 

These  data  thus  provide  the  variation  of  D_  (x)  with  pAx  ;  for 

K 

simplicity  we  will  at  present  place  x^  at  the  origin  (i.e.,  =  Q) 

and  also  assume  that  the  initial  damage  Dq(x)  =  0  .  We  found  that 
-  AV 

D  can  be  approximated  very  nicely  by  means  of  a  least  squares  curve 
R 

fit  as  the  simple  expression 


*  AV,  , 

Dr  (x)  =  a  + 


b  * (l-e~C'pAX) 
pAx 


(63) 


where  a,b',c'  are  positive  constants  for  a  particular  material.  As 
an  example  we  used  the  constant  stress  data  in  [17]  for  AISI  310  stain¬ 
less  steel  at  temperature  T  =  600°C  and  constant  strain  rate  e  = 

o  c 

27.8  x  10  k/sec,  and  obtained  the  values  a  =  0.0015662,  b'  =  0.0013289, 

3  2  2 

and  c'  =  33.0A8296  with  p  =  7.8  Kg/mm  and  A  =  12.6  x  10  mm  . 

From  Equations  (62,63)  we  can  nov  obtain  for  the  combined  material 

parameter  the  simple  exponential  expression 

H(x)  =  ^  (a  +  b  e_CX)  (6A) 

-2 

where  b  =  b'c'  and  c  =  c'pA  in  general,  and  b  =  A. 3918  *  .10 

and  c  =  3.2A80  x  10^  in  the  above  •ample. 
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-  AV 

For  the  above  data.  Figure  3  shows  a  good  agreement  between  D  (>■) 

K 

as  calculated  from  Equation  (63)  and  the  actual  data  points  from  Kef. 
[17].  The  figure  also  contains  a  plot  of  U(x)*Q  in  accordance  with 
Equation  (64).  From  this  curve  [and  Equation  (64)]  we  see  that  H(x)*Q 
reaches  a  peak  of  magnitude  b  superposed  on  top  of  the  constant  value 
a  .  We  may  thus  interpret  the  quantity  b  as  a  measure  of  the  strength 
of  a  local  imperfection.  The  local  inhomogeneous  damage  law  is  now 
given  by  a  combination  of  hquations  (61)  and  (64)  ,  yielding 


1/  ,  ,  -cx,  na+Y  a+<5 

D(x,t)  =  -(a  +  be  )o  t 
U  o 


It  follows  from  Equation  (63)  and  (65)  that  under  the  constant  tensile 
stress  the  local  critical  damage  is 


VO,tCR)  =  DrAV(0)  -  a  +  b  *  DrAV(x) 


where  x  =  x^  =  0  corresponds  to  the  rupture  point.  It  is  clear  from 

Equation  (66)  that  the  critical  value  of  damage  at  the  rupture  point 

is  always  greater  than  or  equal  to  the  average  value.  We  also  have  for 

~  AV 

the  asymptotic  value  of  D 

R 

DrAV)  =  Dr  =  a  (67) 


which  is  in  fact  the  experimental  value  of  the  critical  damage  reported 


in  [1],  From  the  values  of  a  and  b  presented  above  for  AISI  310  stain¬ 
less  steel  at  600°C  and  Equations  (66,67),  we  obtain  for  the  local 

-2 

critical  value  of  damage  D  =  a  +  b  ^  4.5  x  10  and  for  the  experi- 

R 

-3 

mental  critical  value  of  damage  D  =  a  i  1.5  x  10  .  Clearly,  D„ 

R  R 

is  significantly  closer  to  the  theoretical  limiting  value  of  1  than 


0 


As  l)  ■  inis  value  e i  J  tou ic  ne  nr.i>ruvec  luriuer  11  we  nae 
K  t- 

additional  experimental  data  iri  tnc  smaller  interval  U  i  pAx  t  0.1. 

li  we  do  not  place  x^  at  the  origin,  we  can  simply  modify 
Equation  (bij  lor  x  t  x  <  i  as  iollows 


1  .  -c<x-x  .  na+i  u+. 

DU,t;  =  -  ,  La  +  De  o'  ]  j  t 

V  ^ 


Since  Equation  (bfa)  is  valid  onlv  on  one  side  of  x  ,  we  may  extend 

it  further  to  the  case  x  ^  x  bv  assuming  that  tlie  damage  is  symmet- 

o 

rically  distributed  abouc  x^  .  We  thus  obtain  for  U  g  x  (  1 

na+'i  -.it :  , , , 

D(x,t>  =  h  t,x  t  (b9a; 


-C  j  X  -X 


H(x)  =  -  (a  +  be 


For  each  value  of  x  there  is  a  corresponding  rupture  time,  and 
actual  rupture  will  occur  at  tne  minimum  of  these  rupture  times.  The 
rupture  time  for  any  applied  constant  tensile  stress  oq  with  D^fx.i 

assumed  as  zero  may  now  be  obtained  from  Equations  (69)  as 

'  J 

-(na-t-q)  u+ '  /~n\ 

t„  =  -,y-  0  1  '  (/U) 

K  __a+b  o 

As  mentioned  in  Section  jv  ,  .cmpr.’ssion  will  produce  essentially 
no  damage  and  vfe  can  generalize  Equation  (65)  for  both  tensile  and 
compressive  stresses  to  the  form 


IJCX  ,L)  =  -(a  +  be~L  1 "  "  u 


na+Y  a+6  . 

1 

o  o 


If  additional  simplicity  is  recuiie  ,  one  can  approximate  Equation 


(71)  by 


27 


1 

4 


i 


D(x,t)  =  i(a+b[U(x-x  +  y)-U (x-x  -  y)  ]  j-o  ntt+6ttt+6u(o  )  (72) 

0  (.  o  l  o  l  J  o  o 

which  distributes  the  damage  uniformly  over  a  small  region  with  center 
at  xq  and  of  width  d  .  In  this  case  one  obtains  the  same  rupture 
time  as  in  Equation  (70),  but  rupture  will  now  occur  at  any  point  within 
this  region. 

We  shall  now  extend  the  constant  load  solution  lor  homogeneous 
materials  presented  in  [2],  using  the  above  inhomogeneous  formulation. 


4.b  Inhomogeneous  Creep  Rupture  Under  Constant  Load-Large  Deformation 
If  lateral  contraction  due  to  the  Poisson  effect  is  not  ignored, 
the  cross-sectional  area  of  a  bar  under  a  constant  tensile  load  P 

o 

is  clearly  a  decreasing  function  of  time.  As  a  result,  the  stress 
increases  with  time  in  accordance  with 


o(x,t) 


P 


o 

A(x,t) 


(73) 


where  the  area  A(x,t)  may  in  general  also  vary  with  the  axial 
coordinate  x  .  In  order  to  generalize  inhomogeneous  damage  Equation 
(60)  to  the  case  of  variable  stress,  we  first  rewrite  it  in  the  form 


5  j  D  f  x ,  t )  6  y  /  6 

3t  [c(x)e  (x,t)aj  0 

*-  c  -1 


Then  by  analogy  with  the  procedure  employed  for  homogeneous  damage 
(see  12]),  we  postulate  that  for  variable  stress  the  in  Equation 

(74)  may  be  simply  replaced  by  o(x,t).  Integrating  this  result  with 
the  initial  condition  D(x,0)  =  0  ,  we  obtain  the  integral  form 
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t 

-  j  i  ^ 

L»(x,t)  =  Cixjc^x.t)  '  (x,  t  ’) 1  udt 1  l  1,73; 

o 

Equation  (/a)  is  consi.stant  with  the  postulate  that  the  rate  oi  damage 
be. a  function  of  the  creep  strain  rate,  creep  strain,  damage  and 
stress,  i.e. 


~  b(x,t)  =  f".—  <-  c  fx  ,  t  >  ,  >-ctx,t),  iJ(x,t),  o(,x,t)  (76) 

ior  c(x,t)  =  c  Equation  \1ji  learis  reduces  to  Equation  (o‘o)  . 

U 

the  total  lateral  strain  given  as  in  Lourtine  et  al.  [ J6j  by 

•;(.x,tj  =  ~  -b  s  c  (x  ,  t )  (77) 

where  L  is  the  elastic  modulus  and  .  and  v  are  the  elastic 

e  s 

and  steady  creep  Possioa's  ratios  respectively,  lor  simplicity  we 
shall  neglect  transient  creep,  ut  this  can  be  corrected  in  an  approxi¬ 
mate  manner  by  replacing  i’  with  u  ,  a  smaller  fictitious  modulus 
o!  elasticity  (see  .  t  )  .  i  r>  ep  strain  for  an  inhomogeneous  material 

then  1  ollows  t  rom  Equation  'j 8)  !\  as 

l 

•-  c  (x ,  t ;  =  K(x;  .Jtx.t')1'  dt'  (78; 

o 

By  a  combination  of  Equations  t  we  obtain  the  result 


<  L 


(X  ,t; 


vsKtX) 


A(x,t)_  ’ 


t  >  o  (79; 


For  the  strain  rate  we  use  t:,.  logarithmic  icrm  for  finite  deformation 


(x  ,  t ) 


\  xX 


t) 


t»0) 
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where  R(x,t)  is  the  instantaneous  cross-sectional  radius.  This  leads 
to  the  partial  differential  equation  in  A(x,t) 


1  3A/3t  ePo  3A/3t  „„  n ‘'ll’ 

"2  A  ~ E  "a - JSK(X)P0  (Aj 


(81) 


which  at  some  constant  value  of  x  =  may  be  rewritten  in  the 
ordinary  differential  form 

v  p 

vsKJ>ondt  =  An_2  dA  -  2  An_1  dA,  t  >  0  ,  x  =  Xq  (82) 

Integrating  Equation  (82)  with  initial  condition  A(x ,0+)  =  A  ,  we 

o 

obtain  the  solution 


,A 

o 


2nv  c 


co  !  A 


1-1- 


2nv  o 
e  o 


1  -  V  '  'Ac  l~r~  =  1  -  y  ■  y\-,.  -2nv  K  0  t,  t>0 

L  tn-l)E  j  J  (n-l)E  s  o 


(83) 

which  being  valid  for  any  x^  is  in  fact  valid  for  all  x  .  Equation 

(83)  is  a  solution  for  A(x,t)  in  implicit  form,  and  thus  must  in 

general  be  solved  by  numerical  techniques.  However,  we  may  approximate 

——  within  the  brackets  by  unity,  and  obtain  the  approximate  closed  form 
o 

solution 


,/  -  2n(n-l)Ev  K(x) 

A(x ,  t)  j.  _  s  n 

A  ^  I  (n-l)E-2nv  o  o 

o  L  eo 


t  >  0 


(84) 


The  stress  then  follows  from  Equations  (73)  and  (84)  as 


A 


2nfn-l)Ev  Kfx) 


o(x,t)  = 


‘o  A  (x  ,  t ) 


■t  u 


1  _  - 


(n-1) E-2nv 


n 

o  t 
o 


1 

n 


(85) 


Substituting  Equation  (85)  info  Equation  f78)  and  evaluating  the 
integral,  we  also  obtain  the  creep  strain  as 


_>U 


vtiore 


L*. 


;ii-l )  lu-lin*-  i: 

L.‘  O 


1 _ 

iktx)c  11  uc 
o 


°o  ,  which  is  t»iniiiai 


to  iioii'o  ducliie  ru;jtui'i*  criterion.  The  damage  may  now  be  obtained 
from  n(;aat.  icu^  as 


L)^X,  t)  r  C(,X;t  U,t)’ 
c 


on  -  ^ 

:  ':--L - - - •(l-lR(.x)o  Ut)  6n 

°  '^wn  .KU)c  n  _ 
o 


^o7) 


For  the  spec  iai  case  •,  -  in  ..see  discussion  in  Section  j;  ,  Equation 
(,87)  simplifies  to 


D  i  s  ,  t  ,i  C(x)t  U,  t/‘  c  ’  — - — —  ln(l-lK(x)j  n  t) 

c  O'..  n  o 

■  n'.x;c 

a 


l8fc) 


At  each  position  x  and  corresponding  to  a  critical  value  of 

rupture  J  there  is  a  rupture  title  t  (x)  expressed  by  Equation  (o7)  , 

which  is  an  implicit  function  o:  x  and  t]( (x) .  Actual  rupture  will 

occur  at  that  position  in  the  bar  where  there  is  a  minimum  value  oi 

t  (x;  .  Thus,  in  general  A  location  and  time  of  rupture  is  deter- 

mined  by  £  minimi;:  at  id  o;  t.,  :  in  Equation  (67).  In  the  special 

'  and  Atx,  -  K,  tor  which  Equation  (Ebb; 

is  valid,  we  know  a  priori  timt  nature  will  occur  at  x  -  x  ,  and 

o 

that  the  ruptuie  time  l  i.->  tin  loot  o;  the  equation 


.omogeneous  case  D  (x. 


D,  = 


a-t  b 


(t  )  ‘  •  '' 

c  h  o 


Oil 

tl-lRo  U  t.  )  “U  -l1  }  (8b) 

O  k  1 


Wt:  shall  consider  in  tie  in  >: :  ection  inhomogeneous  rupture  under 
constant  I  oati  villi  smaJ  1  del  «.-r •  n. 


4,c  inhomoiceni  ous  Rupture  Under  Constant  Eoad-Sir.il  1  Deformation 


Let  us  dt  l  me  tlu*  area 


,  n  - 
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A(x,t)  =  A  -A(x,t)  ,  A (x , t )  A 

o  o 


(9u; 


where  in  accordance  with  small  deformation  theory  A(x,t)  is  a  small 
quantity.  With  transient  creep  still  excluded,  substituting  Equation 
(90)  into  Equation  (81)  and  neglecting  higher  order  terms  we  obtain  tiie 
linear  partial  differential  equation  with  variable  coefficients. 


( V  c 


eo  l!  5  A/  ot  ,  .  n'A  >  ,  n  n  .  .... 

E~  "  2j~r~  n  VskU',0o  (AJ  VS  k(x)oo  =0,  t>u  (91) 


For  convenience  we  rewrite  Equation  (91)  as 
-,T  .  a  (x) 

a  ~  +  a  (x )A  +  —  =  0  ,  t  >  0  i92) 

lot  l  n 

.  i  A  e  o  l  .  .  ,  ,  ,  n 

where  A  =  —  ,  a^  =  — - j  ,  ana  (x )  =  n  vgh(x)oo  .  At  some 

o 

fixed  position  x  =  x^  ,  Equation  (92)  is  the  ordinary  differential 
equation  with  constant  coefficients 


a  -  2 

a,  -r—  +  A  +  —  =  0,  t  >  0  ,  x=x 
“  n  ’  c 


-A  ' 

1  dt  2 


(93) 


Employing  the  initial  condition  A(x,0+)  =  0  the  solution  is  easily 
obtained  as 


A(x,t)  =  e 


a2 

- t 

al  _1 
n 


t  >  0 


(99) 


which  is  valid  lor  all  x  since  the  result  is  independent  of  the 


choice  of  x 

c 


Substituting  Equation  (94)  into  Equation  (73)  and  neglecting  Lie 


higher  order  terms,  we  obtain  tor  the  stress 

_  a,(x) 

I  3 1  1  i 

o(x,t)=o|l  +  e  -  —  | 

o  |  n  | 

L  -i 


(93) 


The  expression  fur  the  creep  strain  then  :  ollows  iron,  equations  1 7  3 , 9 


ci  S 


a.-,  (x) 


na, 

:c(x,t)  =  R(x).,oU  --^y;l 


(96) 


Finally,  the  damage  For  small  deformation  is  obtained  from  Equations 
(75,95;  as 


,(x) 


v  !  ■  *  ■  7  a  2  V A ) 

D(x,t)  =  E(xjt  (x.t)a  o  {  1-  -f—  t  -t - -  1  -  e 

C  O  Ml;  ta, 


-t  ■ 


( 


(97) 


The  method  discussed  in  Sect  ion  4(b)  for  determining  the  rupture 

time  and  the  location  of  rupture  clearly  also  applies  here.  In  general, 

for  a  given  i)  we  must  minimize  t  (x)  as  prescribed  by  Equation 
K  K 

(97),  tn  order  to  locate  the  position  where  rupture  will  take  place. 

Then  by  substituting  this  particular  value  of  x  back  into  Equation 

(97)  we  may  find  the  value  tT1  .  For  the  special  homogeneous  case 

i.>o(x)  =  0  and  K(.x)  =  K  ,  rupture  will  clearly  occur  at  x  =  xq 

and  the  rupture  time  t,  is  obtai  ted  i rom 

K 


2 

-  -  —  t , 


U„  -  --  ,  uB)  ‘  o  V'(l  -  --  t„) 

R  cx  c  K  o  .  r.  K 

QK 


(1  - 


(98) 


in  the  next  section  we  examine  ct  of  randomness  in  the  material 

parameters  or.  the  rupture  time. 


5.  EFFECT  OF  RANDOM  PARAMETERS  ON  CREEP  RUPTURE 


.  Creep  parameters  are  very  highly  dependent  on  the  microscopic 
imperfection  variations  introduced  in  the  manufacturing  process  and  on 
the  small  temperature  variations  encountered  in  practice,  and  such 
variations  are  very  difficult  to  control.  It  thus  follows  that  this 
lack  of  precise  information  concerning  imperfection  and  temperature 
fluctuations  will  be  reflected  in  considerable  uncertainty  in  predict¬ 
ing  lifetime  to  rupture.  Classical  probability  techniques  were 
employed  by  Cozzarelli  and  Huang  [37-39]  in  the  study  of  steady  creep 
of  structures  with  random  material  parameters.  Huang  and  Valentin  [AO] 
extended  this  work  to  the  study  of  creep  deformation  and  creep  rupture 
in  cylindrical  tubes.  Such  techniques  are  frequently  difficult  to 
apply  unless  some  rather  simple  models  are  introduced  or  some  major 
simplifications  are  made.  In  order  to  reduce  the  difficulty  encountered 
in  such  nonlinear  creep  studies,  researchers  frequently  employ  the 
concept  of  a  reference  stress  at  which  the  strain  rate  is  insensitive  to 
the  creep  power  n  (e.g.,  see  I’ennv  and  Marriot  [41]).  In  a  prob¬ 
abilistic  study  of  creep  rupture  broberg  [33j  observed  that  in  a  In 
versus  In  oo  plot  a  variation  in  n  leads  to  a  pole  at  a  certain 
stress,  i.e.,  at  the  reference  stress  o„  .  Recently  Boyle  [42-44]  has 
also  developed  a  systematic  method  lor  the  application  of  the  reference 
stress  technique  in  the  study  of  creep  rupture  with  random  effects.  in 
this  section  we  will  employ  our  inhomogeneous  strain-dependent  rupture 
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model  to  study  the  <*>  j  eel  01  random  parameters  on  the  rupture  tine  tor 
the  case  of  a  bar  under  constant  tensile  stress  and  a  beam  under  pure 
bending . 


5. a  Random  Inhomogeneous  creep  Rupture  for  Uniaxial  Constant  Stress 
We  consider  first  the  simple  inhomogeneous  constant  stress  case 
which  was  treated  deterministically  in  4(.a).  As  we  previously  noted, 
the  constant  b  in  li(x)  [see  liquation  (59b)]  may  be  considered  a 
measure  of  the  strength  of  a  local  imperfection.  In  this  subsection, 
we  shall  simply  assume  that  the  randomness  in  the  material  occurs 
only  through  random  fluctuations  cf  b  .  Equation  (69)  with  b  a 
random  variable  is  then  written  as 


~ ,  .  ..  i  .  ,  £■  -c  x-x  i,  nu+y  ~a+6 

D(x  ,  t)  =  La  +  b  e  o 1  ]o  t 


(99) 


where  the  carat  denotes  a  random  quantity,  and  it  is  clear  that  now 
the  damage  is  a  random  pmc.-ys.  from  Equation  (99)  we  see  at  once 
that  rupture  will  occur  at  position  x  =  x  and  the  rupture  time  is 
the  random  variable 


t 

k 


1_  _ 1_ 

o  "(nri+Y)7  ’T+6(a+b)~  0+6 


Combining  Equations  (70,100)  we 


o 


fa  -t_C ) 
if 

o 


1 

a+d 


(100) 


(101) 


where  t  and  li  are  the  de 
R  o 

o 


••sir i s‘ 


values  oi  t  and  (a  +  b) 
K 
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for  b  equal  to  its  mean  value  b  .  As  previously  noted,  the  experi¬ 


mental  results  indicate  that  b  is  always  positive,  and  thus  the  random 


variable  b'  (defined  as  b’ 


-)  will  be  distributed  from  zero  to 


infinity.  In  order  to  evaluate  the  effect  of  randomness  in  b'  on  t  , 

K 

we  shall  assume  that  b'  is  a  two  parameters  log  normal  random  variable 
[45],  which  distribution  is  frequently  used  to  describe  randomness  in 
positive  material  constants  (e.g.,  see  [34]). 

The  first  order  density  function  of  £  '  is  accordingly  given  by 

.2 


fc(b') 


(In  b’-n)‘ 
2t^ 


U(b') 


(102) 


/27  t  b ' 

The  two  parameters  t  and  n  in  Equation  (102)  specify  the  mean,  variance 

and  most  probable  value  of  b'  in  accordance  with 

1 


n+ 

E{b '  }  =  n  ,  =  e 
b 


Tn+x*1 


V 


(e  -1) 


>>  =  q— T 


MAX { b  '  }  =  e 

The  first  order  density  function  of 
obtained  as  (see  [46]) 


(103a) 

(103b) 

(103c) 

(Equation  (101))  is  then 


^R 

o  o 


’R 


11  (a+5) 
o 


o  _ _ 

2t^ 


171  [  (HQtR) 


-(a+6) 


-*i(W 


l+a+6 


b'f— ] 

J 


(104) 


and  the  mean  rupture  time  follows  as 
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tR  'i 

Ur  J 


fa+b  b') 
o 


uO  : 


2xi 


db ' 


(105) 


v 2r  tL ' 


which  in  general  must  be  evaluated  numerically. 

By  using  a  Taylor  series  expansion  about  n  ,  and  neglecting 
higher  order  terms,  we  may  obtain  a  simple  closed  form  approximation  for 

.A  i 

Ei - —  >  as 


Ai  ,  J  f 


i .. 


'  I  t„  I  %  H  h  o 
R  J  n  i  i 


r'"r  T1'  a+5 

a+n  e  +b 


.  JL  In 

1  1  1  , ,  n  2 1  "l  a+6  2 

o  J-5 W-  ^  ”a+6_!  (A  e  J 


1  2n+v 


(el  -1) 


(106) 


'  ^R  1 

Another  simple  closed  form  approximation  for  Hi - ?  may  be  obtained 

Ur  J 

through  the  use  of  two-point  est  imates  [ A 7 j  .  By  this  method  the  first 
order  density  function  of  b '  is  approximated  by  two  delta  functions, 
which  have  tiie  same  mean,  second  and  third  central  moments  as  b'  . 
Accordingly,  we  obtain 

ft, 


] 


_  _JL 

,  j  R  [  1  Jr  tl+6  1 

L  2iri[a+bo^'b+Tb;j  +  l^VV*Tb)]  j 


Ur 


(107) 


ft,. 


Table  1  gives  a  comparison  !■'-  — '  as  obtained  by  numerical 

' 

.  o 

integration  of  liquation  (105)  and  by  approximate  formulas  (106,107)  for 
AlSi  310  stainless  steel  at  b00°C,  where  a  =  0.7  and  6  =  0.012,  and 
with  t  =  0.3.  It  shows  that  in  this  case  these  two  closed  form  formu¬ 
las  do  give  a  good  estimate  of  i_(- — j-  .  Also,  Figure  4  gives  curves 
j  L  K  t  *' , 

ot  Li — •  plotted  versus  \  n-i  ,  ’  jlu  stainless  steel  at  600oC 
lCR  1 
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with  r  =  0.5  and  0.8  as  obtained  by  numerical  integration  of  Equation 
(105).  We  see  from  these  curves  that,  for  a  given  value  of  t  ,  values 

f^R  1 

of  n  in  only  a  rather  small  interval  yield  values  of  E-j^— j  in  the 
vicinity  of  1.  In  the  next  subsection  we  treat  a  more  complex  example, 
and  also  allow  several  parameters  to  be  random  variables. 


5.b  Random  Creep  Rupture  in  a  Beam  Under  Pure  Bending 

Consider  the  same  problem  treated  deterministically  in  3(b),  but 

now  with  a  random  inhomogeneous  damage  coefficient  and  random  creep 

parameters.  For  the  reasons  previously  discussed,  we  shall  first  rewrite 

Norton's  law  in  terms  of  a  reference  stress  o  in  the  form 

K 

ec  =  eR|-|n  sgn(a)  (108) 


•  A.  • 

where  e  and  n  are  in  general  random  variables,  and  where  e  is 
R  K 

the  strain  rate  at  the  reference  stress  .  Note  that  at  o  =  o_  the 
strain  rate  e  is  unaffected  by  random  fluctuation  in  n  (see  Refs. 
[37]  for  a  detailed  discussion  of  Equation  (108)).  Proceeding  as  in 
the  deterministic  case  we  obtain  governing  equations 


y 


K 


M  =  o, 


n  3  sgn(<) 


J  = 


2n 

2n+l 


wh2hn 


1 

„  M  |  y  |  ^  x 

Q  =  ‘‘■i-  sgn(y) 

j 

=  2n+l  M_ 
°max  2n  wh2 


(109a) 

(109b) 

(109c) 

(109u) 

(109e) 


3b 


Also,  the  damage  law  wrtn  D  (>■;  =  0  will  oe  taken  as 

o 

D(x,t)  =  C(x)  tV.t)  c  yc6  (110) 

IUdA 

where  C(x)  is  a  random  process  in  x  . 

We  shall  assume  that  C(x)  takes  the  similar  form  as  Equation 
(6yb)  ,  i.e., 

.  -  -C1  x_x 

£(x)  =  a  +  b  e  1  o'  (111) 


but  now  both  a  and  £  may  be  random  variables  due  to  the  material 
imperfections.  Then  the  damage  is  a  random  process,  expressed  in  terms 
of  four  random  variable  material  parameters,  obtained  by  substituting 
Equations  (108,111)  into  Equation  (110),  giving 

o_.  j  cxn+y 


D(x,t)  =  (£  -r  fi  e-c'x-V):..  V1  — '  t 

K  R  '  °R 


a+6 


(112) 

Although  there  are  several  expressions  for  the  reference  stress  o 

R 


given  in  [43],  they  differ  very  little  from  the  simplest  expression 
M 

°K  =  wh2  f°r  tlie  interval  3  $  n  S  11  .  Therefore,  for  convenience 
we  choose  Cj,  =  ~y-  in  Equation  (112)  and  obtain  a  simpler  form  for 
the  damage  law  as 


„  -c : x-x 

0(x.t)=(a+Ce  1  o'). 


y  2n-i  1 1  ur.+v  ^a+6 

— :  ■—  t 


(113) 


1' .  om  Equation  (113;,  it  is  clear  that  rupture  will  occur  at  x  =  xq 

(if  Dq(x)  =  3;,  and  the  initial  rupture  t irae  (i.e.,  the  time  at  which 

the  beam  begins  to  rupture)  is  it.  random  variable 

J 


II, 


ill  -  .  ill  tv ) 


-  . ,  + 1 


’HI 


Ha  3  6>;R%Rn  -  3 


(114) 


Assuming  that  a  , 


and  n  are  independent  random  variables. 
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we  may  find  an  expression  for  the  mean  of  the  initial  rupture  time  as 


R1 

tRI  f  Bo 

Q  —00  —00 


(iTb)a+"  faU)  fCCb)  dadb 


(1  a+6  f-  (-  •  (££-) 

c  E  K  lx  J  2n 


an+y 
2n+l.  a+6 


f  * (n) dn  (115) 

n 


.  an  +v 

1  a  0  _  o  ’ 

i  „  ,  1  .  a+6,1  ,a+5.^no+  .  a+6  ,  .  , 

where  B  =  ( — ~r— )  (y— )  (~~ - )  and  t„T  is  the  determin- 

o  a  +b  e„  2n  R1 

o  o  R  o  o 

istic  value  of  £  when  a,  b,  £  and  n  are  equal  to  their  mean 

hi  K 


values  a  ,  b  .  e„  and  n 
o  *  o  R  o 

o 

As  in  last  section,  we  may  approximate  the  mean  of 
Taylor  series  expansion,  yielding 


~RI 

R1 


by  a 


Jk_l iJ  ( 

Vri 

o 


_  JL.  „  -  (  1  +2 ) 

s  a+.'  1  1  ,  1  |  /  ,  N  a+6  ,  i  , 

n  +n,  )  +  7~rr\7  +  ~rrl  (r',+r‘,J  (i~  + 

a  b  2  (u+6)z  a+6,  a  b  a 


2>  1  «  J  ,  \  a+6  if/  u  ■, 2  ,  a  ].  .  'a+6 

vf  1(%t>  +ij_W  +  r+sj > 


“V1 
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2  4-1  —  1  T 

n  ^  a+6  1_  d^  I  , 2n+l.  a+6 

2n_  j  +  2  dn2 2  2n  ' 


1  n=  r  n 
1  n 


(lib) 


We  may  also  use  the  two-point  estimate  technique  to  obtain  an  alterna¬ 


tive  simple  approximate  expression  ior  E 


.j^KI  j 


IsT/ 


•R1 


E^U1  1 


F  Al[(VTa)  +  (r'b+1b)]  a+^(na+Ta)  +  (VTb)] 


1 

a+6 


+  [  (na-Ta)+(nb+Tb) )  a+6+[  (na-ia)+(r,b-ib)  ] 


1 

a+6  ^ 
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b.  BRIEF  SUMMARY 

The  continuum  mechanical  and  the  physical  metallurgical  approaches 
employed  in  various  previous  deterministic,  homogeneous,  strain- 
independent,  uniaxial  creep  rupture  theories  were  first  discussed 
in  some  detail.  Previous  multiaxial  creep  rupture  theories  were 
also  briefly  discussed,  with  the  damage  treated  either  as  a  scalar  or 
as  a  tensor  quantity.  Then  the  deterministic,  homogeneous,  uniaxial 
tensile,  strain-dependent  creep  damage  theory  developed  by  Belloni, 
Bernasconi,  Cozzarelli  and  Piatti  from  both  the  continuum  mechanical 
and  physical  metallurgical  approaches  was  studied  in  considerable 
detail. 

Based  on  previous  experimental  results  obtained  at  J.R.C.  Ispra 
it  was  postulated  that  no  damage  is  produced  by  compression,  and 
using  this  postulation  the  Belloni  et  al.  strain-dependent  creep 
damage  theory  was  generalized  to  include  compressive  stress.  A 
multiaxial  strain-dependent  creep  damage  theory  was  also  developed 
by  a  simple  extension  of  the  uniaxial  theory  with  the  damage  treated 
as  a  scalar.  The  problem  of  a  beam  under  pure  bending  was  then 
studied  using  this  strain-dependent  creep  damage  theory.  At  the 
instant  of  initial  rupture  a  failure  front  is  formed  at  the  top  out¬ 
side  fiber  of  the  beam,  and  then  this  failure  front  propagates 
through  the  entire  cross-section,  resulting  in  a  complete  rupture  of 
the  beam. 

An  analysis  of  creep  rupture  data  from  J.R.C.  Ispra  for  speci¬ 
mens  of  various  lengths  at  constant  temperature  and  under  constant 
tensile  stress  resulted  in  the  postulation  of  a  local  inhomogeneous 
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ext  ending  previous  results  as.  ..  or.  homogeneous  ti.ecry. 

Also,  t  :.e  jnhonoceuet  us  .->t»  a  i:. -dependent,  creep  rupture  model 
was  used  to  examine  the  >•:  feet  random  put  arte  t  e  rs  on  rupture  i  ice 
using  classical  probub  ;  1  i  c>  tec;  :.  J  ..-ties.  it:  cite  case  oi  a  bar  under 
constant  tensile  stress,  the  nia.  :  itude  of  ti.e  peak  ot  the  combined 
material  parameter  va-  t;»ui.-d  a-  a  randoi;  variable.  Vhe  mean  value 
of  the  rupture  i  into  was  then  obr.his  1  f ror.  numerical  integration  anu 
from  two  simple  closed  form  at.  ox  in-.at  ions  which  were  in  g-oc  agree¬ 
ment  with  the  results  oi  the  nun.  i  ical  integration.  This  ir.eau  value 
demonstrated  tii.it  the  rupture  t  •  • is  highlv  cepeiident  on  random 
imperfection  variations.  Ki  a!  ,  .  beau  under  pure  bending  witli 
several  parameters  chosen  as  random  variables  was  reexamined  using 


a  similar  approach 
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Figure  3,  Average  Damage  at  Rupture  D^v  and  Combined 
Material  Parameter  H(x)-Q  vs.  Mass  p Ax 
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